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Abstract 

The production of the hghtest Higgs boson, h^, in association with the hghtest stop, 
ti, in the Minimal Supersymmetric Standard Model, is an interesting channel, as its 
cross section can be higher than Higgs production in association with top quarks. Fur- 
thermore, the iii\h^ production rate is highly dependent on various Supersymmetric 
parameters. The mass, mixings, couplings and production cross sections relevant to 
the tit*h^ channel have been studied in the past. Here, we complement these anal- 
yses by performing a thorough decay study. We conclude that there is some scope 
for extracting this channel at the Large Hadron Collider, for suitable combinations of 
the Supersymmetric parameters, in several different decay channels of both the stop 
quarks and Higgs boson, the most numerically promising being the signature involving 
charm decays of the stops and bottom decay of the Higgs. This, in fact, remains true 
in the Supergravity inspired Minimal Supersymmetric Standard Model. 
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Introduction and motivations 



The Standard Model (SM) of particle physics has been extremely successful at explaining 
experimental data to date. It has, however, failed to provide answers to fundamental prob- 
lems, such as the unification of the gauge couplings at high energies, as well as the difference 
between such a unification scale and that of Electro- Weak (EW) interactions (known as the 
hierarchy problem). Introducing a new fundamental symmetry of Nature between bosons 
and fermions, known as Supersymmetry (SUSY), solves both these problems. If SUSY ex- 
ists at the EW scale, it should be discovered at the Large Hadron Collider (LHC). In the 
Minimal Supersymmetric Standard Model (MSSM), all known particles will have Super- 
symmetric partners, known as sparticles. Additionally, in the MSSM, two complex Higgs 
doublets are required to avoid anomalies, which in turn have Supersymmetric partner states. 
EW Symmetry Breaking (EWSB) results in five Higgs bosons in the MSSM. The lightest 
neutral Higgs state of the MSSM, h^, is likely to have a mass of less than 135 GeV, so it 
is well within the kinematic reach of the LHC. Current experimental limits on the mass of 
the sparticles that could enter decay channels are such that (with the possible excep- 
tion of invisible decays, e.g., into two neutralinos, the SUSY partners of the SM neutral 
gauge/Higgs bosons, which are however small [HEIE]) the decay modes that are interest- 
ing phenomenologically are those involving SM particles (henceforth we will refer to these 
as SM-like channels, likewise for the production modes involving couplings to SM particles 
only). However, it does not follow that decay modes and possible discovery channels 
should be the same as for a low mass SM Higgs state (hereafter denoted by H). In fact, 
Supersymmetric effects can affect both the Higgs production and decay rates. 

The tit\h^ production channel is interesting for two main reasons: 1. it provides a 
production mechanism for the lightest Supersymmetric Higgs boson, in addition or as an 
alternative to the SM-like channels; 2. its production rates are strongly dependent on the 
MSSM parameters. For M/jO < 135 GeV, the dominant decay mode is — > 66, however, at 
the LHC the inclusive bb decay is not visible due to the large QCD background. The 

— > 77 decay mode is the most probable LHC discovery channel for a Higgs boson at these 
masses [1]. The — 77 decay, however, can be affected by Supersymmetric effects [5], in 
such a way that cancellations between the SM particle loops and those involving their SUSY 
partners can occur over large regions of the MSSM parameter space. This may deplete the 
liP decay rates below experimental observabilit}@. This typically occurs for light stop 

quarks, with a mass of order of that of the top quark or below. In turn, this means that 
one could re-establish the lost SM-like di-photon signal by exploiting Higgs production in 
association with (light) stop quarks. Even when associated Higgs production and decay has 
a lower cross section than the inclusive di-photon channel, the extra final state particles can 
render this SUSY channel visible over the background by providing useful triggers, thereby 
possibly also enabling — > bb detection. 

^This applies not only to the exclusive ^"fi^vi channel but even more so to the inclusive mode, as a 
similar mechanism onsets in the suppression of the gg production channel. 
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Under the above circumstances, one can elaborate on point 1 as follows. Take for example 
the SM counterpart of iit\hP , i.e., Higgs boson production in association with top quarks, 
tth. This channel has been studied by both ATLAS [6] and CMS |7]. Here, the associated 
top quarks decay to 6-quarks and one of the t-quarks decays leptonically and the other 
hadronically. Since the dominant decay of the H state is ^ — *^ bh, a four 6-quark final state 
is produced. This is much more visible over the QCD background than the inclusive two 
6-quark final state. The scope of this channel in the SM is now questionable though, as 
its feasibility is not certain [7]. However, in the MSSM, when there is a large stop mixing, 
the coupling of the lightest stop to the lightest Higgs boson can be larger than the top 
Yukawa coupling. This means that in the MSSM the titlh^ cross section can be larger 
than the tih one. In addition, for stop masses suitably larger than mj, the ti txi decay 
could occur with a Branching Ratio (BR) nearing one. The final outcome would then be 
a signature similar to that emerging in the SM from tth, possibly with larger rates and an 
additional handle for background rejection, the larger missing transverse energy (due to the 
two neutralinos escaping detection) with respect to the SM case (when it is only due to one 
neutrino produced in the leptonic top decay). The tit\h^ channel therefore has the potential 
to be an interesting channel for Higgs discovery in the MSSM and this is the reason why it 
has received considerable attention in the literature [8],[9]. 

The iiilh^ channel can also provide information on the underlying Supersymmetric sce- 
nario, point 2 above. If SUSY exists at the EW scale, particularly in the form of light 
squarks, it should be discovered at the LHC in direct tiil production (i.e., without subse- 
quent emission). However, information on the particular Supersymmetric model and its 
parameters will be hard to determine from this process (primarily because it proceeds via 
QCD interactions, which are not modified with respect to the SM by SUSY, so that the only 
SUSY parameter measurable is here mi^). The tit\h^ cross section, however, is dependent on 
several low energy MSSM parameters. These enter the tit\h^ cross section through both the 
lightest stop mass and its couplings to the lightest Higgs boson. Measuring the tit\h^ cross 
section could therefore provide us with decisive information on the SUSY-breaking scenario 
and its parameters. Importantly, such low energy SUSY parameters can directly be related 
to those defining the Supergravity inspired MSSM (MSUGRA), which is one of the simplest 
possible scenarios in the choice of the soft SUSY-breaking parameters at the Grand Unifica- 
tion (GUT) scale. In MSUGRA, the scalar masses, scalar couplings and gaugino masses are 
all equal (separately) at the GUT scale. The remaining two parameters needed to describe 
soft SUSY breaking in MSUGRA are the soft Higgs mixing parameter, 6, and the Higgs 
mixing parameter of the two doublets, /i. However, using constraints from the known value 
of the Z mass, these two parameters can be calculated from two EW scale parameters: tan /5 
(the ratio of the Vacuum Expectation Values (VEVs) of the two Higgs doublet fields) and 
the sign of /i. This means that the over one hundred MSSM parameters can be reduced to 
the following five MSUGRA ones. 

1. Mq, the universal scalar mass, defined at the GUT scale. 
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2. Ml, the universal gaugino mass, defined at the GUT scale. 
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3. Aq, the universal trilinear coupling, defined at the GUT scale. 

4. tan P, defined at the EW scale. 

5. sign(/i), the sign of the Higgs mass parameter, defined at the EW scale. 

The tit\h^ cross section depends at lowest order on tan/3 and sign(/i) defined at the EW 
scale and - through the Renormalisation Group Evolution (RGE) equations of the model - 
also on the Mq, Mi and Aq GUT parameters. 

In this paper we examine possible decay channels and test their suitability for detection 
at the LHC. In the next section, we review existing literature and update old results by 
using more recent simulation software. We also perform new scans over the MSSM and 
MSUGRA parameter spaces, where both production and decay rates are determined. Finally, 
we conclude in the last section. 

The iii\h^ rates in the MSUGRA parameter space 

Since the publication of [H [9], the experimental constraints have been revised, and the 
theoretical models have been improved. In this, and the following, section we describe the 
updates to these results, using the latest simulation packages, primarily those that determine 
the MSSM and MSUGRA spectra of masses, couplings and BRs. Here we use ISASUSY 
and ISASUGRA [10] and both were interfaced to the Supersymmetric version [11] of the 
HERWIG Monte Carlo (MC) program [12], where the production process (tested against 
FeynHiggs [13]) is implemented alongside the (matrix element modelled) decays of the final 
state objects. 

There are previous papers that describe studies of Higgs production in association with 
light squark pairs. The first is Ref. |8j. Here, scans over the five MSUGRA parameters were 
carried out (for all processes of the type squark-squark-Higgs, though in this paper we only 
concentrate on the iii\h^ case). The parameter space explored was defined as follows: 

• Mo = Ml =150 GeV. (These low stop and gaugino masses were chosen to give low stop 
and Higgs masses in order to enhance the cross section.) 

• Aq was set to —300, 0, 300 GeV, in order to investigate its behaviour as it changes sign. 

• Both values of sign(/i) were investigated. 

• tan/5 was scanned over values between and 40. 

Figure [1] shows the tit\h^ cross section results from the scans over this MSUGRA space. 
These cross sections are sizably different from those presented in the original paper. Different 
codes were used, however, we have traced back this difference to an error in the convention 
used by [H] in the definition of the stop mixing angle, 9t- While this error is of no consequence 
for the calculation of the Higgs (the angle enters here through the square of its sin function) 
and stop (which is independent by definition) masses, it can have a sizable impact in the 
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Figure 1: tit^h^ cross section values resulting from the variation of MSUGRA parameters. 



tit\h^ coupling, as the latter can be written as 



9htiho = cos 2P ( - cos^ 0t - - sin^ 9w 
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where My^ijz is the mass of the WjZ gauge boson and Ow is the weak mixing angle. 

The cross sections given in figure [T] are generally lower than those presented in Ref. [8]. 
Figure [T] shows that, in general, the cross section decreases with increasing tan/3. Further- 
more, it is highest for Aq = —300 GeV and /i > 0. This is to be expected, as it will give the 
highest value for \At\, as the negative gauge corrections increase the magnitude of a negative 
Aq. There are also minima in the cross section for two of the scans, which were not seen in 
the scans presented in [S] , and are a result of the new calculation of the stop mixing angle. 

In [S], two scans over Mi were also investigated. The two scans are defined in table [Tj 
as the "general low mass spectrum" and the "general high mass spectrum" . 





Mo (GeV) 


Ml (GeV) 


Ao (GeV) 


tan P 


sign// 


Low mass spectrum 


200 


100 - 350 


-300 


35 


-1 


High mass spectrum 


500 


100 - 500 


-900 


35 


-1 



Table 1: The MSUGRA parameters corresponding to the high mass and low mass spectra. 



Figure 2(a) presents the MSUGRA scan for the low mass spectrum whilst figure 2(b) 
that for the high mass one. Figure [2] shows that, in both scans, the cross sections decrease 
as Ml increases. The RGE equations that give the values of the stop and Higgs masses at 
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the EW scale include correction terms from gauge interactions. The latter are dependent on 
Ml, giving an increase in both the Higgs and stop masses as Mi gets larger. These higher 
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(a) titlh^ cross sections for a low mass MSUGRA spec- (b) tiilh'^ cross sections for a high mass MSUGRA spec 
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Figure 2: cross section values resulting from the variation of Mi . 



masses lead to a phase space suppression of tit\h^ production. The cross sections calculated 
here lead to the conclusion that it is unlikely that the tit\h^ channel will be visible in the 
MSUGRA regions that were considered in [8]. 

However, the region of MSUGRA parameter space investigated in Ref. [5] is more promis- 
ing. This is characterised, with respect to the region investigated in [S], by a very large 
negative value of Aq. The region investigated is defined in table [2l The range of Mq values 
used was chosen to give a range of m^-^ values that are likely to be detectable at the LHC. 



Mo 


(GeV) 


Ml (GeV) 


(GeV) 


tan (3 


SigUyU 


700 


- 1200 


300 


-2000 


2.5 and 30 


1 



Table 2: The region of MSUGRA space investigated in Ref. [9]. 



Figure [3] shows how the tii^h^ cross section varies with mi^ over the above area of the 
MSUGRA parameter space. The cross section results in figure [3] show rather good agreement 
with those presented in [9|. As expected, the cross section decreases with decreasing m^^. 
Nonetheless, here, the production rates are high enough, so that, for a large range of mi^ 
values, they could be visible at the LHC. The production rates are not significantly affected 
by the sign of /i. 

The final state occurring after the decay of the stop quarks and Higgs boson entering 
these events also needs to be considered, before drawing any firm conclusions. The most 
common decay modes for the stop are to a top quark and a neutrahno (as mentioned already), 
a bottom quark and a chargino or a charm quark and a neutralino. Figure H] shows how 
the decay modes of the stop change over the range of m^-^ considered in figure [31 Typical 
masses for the two lightest neutralinos in this region are 130 GeV and 250 GeV, with the 
lightest chargino also around 250 GeV. At m^-^ > 300 GeV the stop is most likely to decay 
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Figure 3: tit\h^ cross section values resulting from the variation of m^^, over the parameter 
space defined in table [2l 




(a) ii BRs for tan/? = 2.5. 
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(b) ii BRs for tan/3 = 30. 



Figure 4: ti BRs resulting from the variation of m^-^, over the parameter space defined in 
table El 
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Figure 5: o'tj*ho for the common stop decay modes, over the parameter space defined in 
table [3 

to a top quark along with the lightest neutralino. For m^^ < 250 GeV the most likely decay 
is to the charm quark and a neutralino. A decay to the bottom quark and lightest chargino 
is most likely for 250 GeV < m^-^ < 300 GeV . There is also a small area of this parameter 
space, for m^^ ~ 250 GeV, where the stop is most likely to decay to a charm and the second 
heaviest neutralino, due to a combination of coupling and phase space effects. For m^-^ < 130 
GeV, the stop is the Lightest Supersymmetric Particle (LSP) and therefore will not decay 
(as long as -R-parity is conserved). The value of tan/5 has a rather negligible effect on the 
stop decays. 

Figure [5] shows the variation in the cross sections for tit\h^ , over the parameter space 
defined in table [21 for the two stops decaying in each of the common stop decay modes. 
It can be seen from figure E] that the highest cross section for tit\h^ occurs when the stop 
decays to a charm quark and a neutralino. If the tit\hP channel can at all be accessed at the 
LHC, then it is likely that this will occur in the parameter region where the stop decays to 
a charm quark and a neutralino. 

The most common decay modes for the Higgs boson in the area of parameter space 
considered in [S] are those into bottom or charm quarks or to taus. The BRs to these decay 
modes do not vary significantly over the parameter space considered. This is because the 
Higgs mass only varies between approximately 97 and 120 GeV, much less than the variation 
in the stop mass. The decay fraction of the Higgs to bottom quarks is around 85%, to taus 
around 5% and to charm quarks around 3%. The high — > hb decay fraction, combined 
with a 6-tagging efficiency of around 50%, may enable one to extract decays produced at 
the LHC in association with stop-antistop pairs. 

Having established that the cross sections and BRs are such that it could be possible 
to access the tit\h^ channel at the LHC in the MSUGRA scenario of Ref. [9J, we study an 
expanded region of parameter space. To this end, additional scans were carried out using 
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Figure 6: tit\h^ cross section values resulting from the variation of MSUGRA parameters. 

the adaptive scan program adScan [H]. The adScan program uses adaptive integration 
algorithms to identify regions of interest in a multidimensional scan. It can therefore scan 
over a wide range, and large number, of parameters faster than a uniform scan. The values 
of the MSUGRA parameters were varied as shown in table [31 These values were chosen 
to provide a wide Supersymmetric spectrum, but one that is theoretically viable (i.e., only 
non-tachyonic particles are present, EWSB is achieved, etc.). Although sign(/i) was initially 
kept constant and equal to 1, further investigation showed that the value of sign(/i) has little 
effect on tit\h^ production. 



Mo (GeV) 


Ml (GeV) 


(GeV) 


tan f3 


sign(/x) 


- 1000 


- 1000 


-2500 - 2500 


2 - 30 


1 



Table 3: The region of MSUGRA space investigated using the adScan program. 

Figure [6] shows how the tit\h^ cross section varies for this multi-dimensional scan over 
the MSUGRA space. In particular, it shows that the cross sections for tit\h^ with ti —>■ cXi 
can be as high as 1 pb. For the ti bxt decay channel the cross sections reach a maximum 
of around 100 fb while for ti — > txi the maximal value is only around 10 fb. 

Production cross sections for the MSUGRA tit\h^ signal at the LHC are therefore highest 
for ti cXi (and charge conjugate) and bb. 

The iit\h^ rates in the MSSM parameter space 

While MSUGRA is an appealing representation for a SUSY-breaking scenario, it need not 
be the one realised in Nature. We have therefore scanned over a generic (non-universal) 
MSSM parameter space, defined at the EW scale, where the important parameters for stop- 
stop-Higgs hadro-production are the following. 
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Figure 7: titlh cross section values resulting from the variation of MSSM parameters. 

1. tan/3, the usual ratio of the VEVs of the Higgs fields. 

2. fi, the higgsino mass parameter. 

3. At, the stop trilinear coupling. 

4. m^^, the left-handed stop mass. 

5. rni^, the right-handed stop mass. 

6. M40, the CP-odd Higgs mass. 

As the values of these parameters are limited in the MSUGRA scenario, it may be possible 
that in the general MSSM there are areas of parameter space where the tit\h^ cross section 
is higher than in MSUGRA. 

Multi-dimensional scans over the MSSM parameter space were again carried out using 
adScan. The values of the MSSM parameters were varied as shown in table HI Again, as 
for the MSUGRA scans, values were chosen to provide a wide, but theoretically viable. 
Supersymmetric spectrum. 



mi, (GeV) 


mn (GeV) 


At (GeV) 


tan/9 


/i (GeV) 


niA (GeV) 


- 1000 


- 1000 


-2500 - 2500 


2 - 30 


-2500 - 2500 


- 1000 



Table 4: The region of MSSM space investigated using the adScan program. 



Figure [7| shows how the tit\h^ cross section varies for this multi-dimensional scan over 
the MSSM space. In particular, it shows that the cross sections for tit\h^ with ii — > ex? can 
be as high as 100 pb. For the ti bxt decay channel, the cross sections reach a maximum 
of around 10 pb and for ti txi the maximum cross section is around 100 fb. We therefore 
conclude that within this region of the MSSM parameter space, as for the MSUGRA study, 
charm quark decays of stop quarks provide the highest production cross sections for the 
iiilh° channel at the LHC. 

The MSUGRA and MSSM parameter scans above have not enforced experimental limits 
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Figure 8: tit\h^ cross section values resulting from the variation of MSSM parameters after 
the experimental constraints in table have been taken into account. 



on the (s)particle masses. Table [5] shows these for the sparticles contained in the possible 
tit\h^ decay chains. All these limits are taken from LEP2 and Tevatron searches. A summary 



Sparticle 




h 


xt 


X? 


Experimental Limit (GeV) 


93 


95 


103 


39 



Table 5: Experimental limits on the masses of the sparticles relevant to this study. 

of the experimental limits on (s)particle masses can be found in [T5] . 

In order for the experimental limits to be imposed in a comprehensive way, all sparticle 
mass limits should be considered, as each parameter point in the MSSM or MSUGRA space 
can yield SUSY states not relevant to the tit\h^ cross section yet being inconsistent with 
experiment. Clearly, it is possible to alter the masses of (most) other particles without 
changing the tit\h^ cross section. This can be done by altering the values of any of the 
MSSM parameters, bar the six that affect the tit\h^ channel. 

Performing a scan respecting the experimental limits given in table [5] results in the tit\h^ 
cross sections shown in figure [HI This shows that the cross sections for tit\h^ with ti cx\ 
or ti — * hxi can reach a few pb. For the ti txi decay channel the cross sections reach a 
maximum of around 50 fb. 



Outlook and summary 

In summary, mult i- dimensional scans were carried out over the MSUGRA and MSSM pa- 
rameter spaces. It was found that the highest iiilh^ cross sections correspond to low stop 
masses. In these areas the stop is most likely to decay to a charm quark and the lightest 
neutralino. Considering that the /i° almost always yields bb pairs, the complete final state in 
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presence of these decays will be two charm quarks, two bottom quarks and missing energy. 
The final state for the areas where the stop decays to a bottom quark will be four bottom 
quarks, missing energy, plus two pairs of jets and/or leptons from the chargino decays. The 
same final state will occur for the stop to top quark decay, with the additional jets or leptons 
originating from the W boson produced in the top decay. There is also the possibility of 
mixed stop decays, leading to even more variation in the final states. The diversity in the 
final states means that different triggers, backgrounds and selection algorithms will need to 
be considered for each decay channel. Further detailed studies will therefore be required in 
order to ascertain the visibility of the tit\h^ channel at the LHC 
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